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Abstract 

An approach to analyse the properties of a particle system is to compare it with 
different processes to understand when one of them is larger than other ones. The 
main technique for that is coupling, which may not be easy to construct. 
We give a characterization of stochastic order between different interacting particle 
systems in a large class of processes with births, deaths and jumps of many parti- 
cles per time depending on the configuration in a general way: it consists in checking 
inequalities involving the transition rates. We construct explicitly the coupling that 
characterizes the stochastic order. As a corollary we get necessary and sufficient con- 
ditions for attractiveness. As an application, we first give the conditions on examples 
including reaction-diffusion processes, multitype contact process and conservative dy- 
namics and then we improve an ergodicity result for an epidemic model. 
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1 Introduction 



The use of interacting particle systems to study biological models is becoming more and 
more fruitful. In many biological applications a particle represents an individual from a 
species which interacts with others in many different ways; the empty configuration is 
often an absorbing state and corresponds to the extinction of that species. An important 
problem is to find conditions which give either the survival of the species, or the almost 
sure extinction. When the population in a system is always larger (or smaller) than the 
number of individuals of another one there is a stochastic order between the two processes 
and one can get information on the larger population starting from the smaller one and 
vice-versa. 

Attractiveness is a property concerning the distribution at time t of two processes with 
the same generator: if a process is attractive the stochastic order between two processes 
starting from different configurations is preserved in the time evolution (see Section [2. ip . 

The main technique to check if there is stochastic order between two systems is cou- 
pling: if the transitions are intricate an increasing coupling may be hard to construct. 
The main result of the paper (Theorem 12. 4^ Section 12. ip gives a characterization of the 
stochastic order (resp. attractiveness) in a large class of interacting particle systems: in 
order to verify if two particle systems are stochastically ordered (resp. one particle system 
is attractive), we are reduced to check inequalities involving the transition rates. 

A first motivation is a general understanding of the ordering conditions between two 
processes. The analysis of interacting particle systems began with Spin Systems, that are 
processes with state space {0, 1}^ . We refer to [llj and [12j for construction and main re- 
sults. The most famous examples are Ising model, contact process and voter model. These 
processes have been largely investigated, in particular their attractiveness (see [HI Chap- 
ter III, Section 2]). Many other models taking place on X'^ , where X = {0, 1, . . . M} C N, 
that is with more than one particle per site, have been studied. Reaction-diffusion pro- 
cesses, for example, are processes with state space (hence non compact), used to model 
chemical reactions. We refer to [H Chapter 13] for a general introduction and construction. 
In such particle systems a birth, death or jump of at most one particle per time is allowed. 
But sometimes the model requires births or deaths of more than one particle per time. 
This is the case of biological systems with mass extinction ( [l7j , [18] ) , or multitype contact 
process ([6], [7j, [H], [K]). A partial understanding of attractiveness properties can be 
found in [20] for the multitype contact process. A system with jumps of many particles 
per time has been investigated in [HI Theorem 2.21], where the authors found necessary 
and sufficient conditions for attractiveness for a conservative particle system with multiple 
jumps on with misanthrope type rates. 

Those examples and the need for more realistic models for metapopulation dynamics 
systems (|9]) has led us to consider systems ruled by births, deaths and migrations of more 
than one individual per time with general transition rates to get an exhaustive analysis 
of the stochastic order behaviour and attractiveness. Our method relies on [8], that it 
generalizes. 

The main applications allow to investigate the ergodic properties of a process. A pro- 
cess is ergodic if there exists a unique invariant measure to which it converges starting from 
any initial configuration: if the process is attractive, it is enough to check the convergence 
of the largest and the smallest initial configurations. This is a first application of Theorem 
12.41 In Section 12.2.51 we combine attractiveness and a technique called criterion (see 
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[3]) to get ergodicity conditions on a model for spread of epidemics, either if there is a 
trivial invariant measure or not. 

For many biological models the empty configuration is an absorbing state and the 
main question is if the particle system may survive, that is if there is a positive probabil- 
ity that the process does not converge to the Dirac measure Sq concentrated on 0, which 
is a trivial invariant measure. In order to prove that a metapopulation dynamics model 
(see [3]) survives, we largely use comparison (therefore stochastic order) with auxiliary 
processes: this is a second application of the result. Instead of constructing a different 
coupling for each comparison, we just check that inequalities of Theorem 12.41 are satisfied 
on the transition rates. Moreover the main technique we use to get survival is a compar- 
ison with oriented percolation (see [Bj), and attractiveness is a key tool in many steps of 
the proofs. 

The survival of a process does not imply the existence of a non trivial invariant mea- 
sure: one can have the presence of particles in the system for all times but no invariant 
measures. If the process is attractive and the state space is compact, a standard approach 
allows to construct such a measure starting from the largest initial configuration: this is 
the third application. Once we get survival, we use this argument to construct non trivial 
invariant measures for metapopulation dynamics models. In Section 12.2.41 we introduce 
a metapopulation dynamics model with mass migration and Allee effect investigated in [3]. 

The transition rates of the particle systems we analyse in this paper depend on two 
sites X, y, on the number of particles at x and y and on the number of particles k involved 
in a transition: they are of the form b{k,a, f3)p{x,y), where a and f3 are respectively the 
number of particles on x and y and p{x, y) is a probability distribution on TJ^ given by a 
bistochastic matrix (we require neither symmetry nor translation invariance). Moreover 
we allow birth and death rates on a site x depending only on the configuration state on x. 

In other words we work with three different types of transition rates: given a con- 
figuration T], on each site y we can have a hirth {death) of k individuals depending on 
the configuration state rj{y) on the same site y with rate -P^^^) ^^niy)^ ^"^^ depending 
also on the number of particles on the other sites x ^ y with rate r){y)P^^^ 
iYlx ■^^(y)r]{x)P^y^^^^' consider a death rate R~^yy^f^^-^p{y , x) instead of the more nat- 
ural R~^^^ ^i^^^p{x , y) to simplify the proofs and since we are interested in applications 
given by a symmetric probability distribution p(-,-). This represents a possible different 
interaction rule between individuals of the same population and individuals from different 
populations. We can have a jump of k particles from x to y with rate T^j-^^ ^^^^p(x, y), 
which represents a migration of a flock of individuals (see Section 12. ip . We require that 
the birth/death and jump rates differ only on the term 6(A;, a, that is the conservative 
and non-conservative rates depend on the same probability distribution p{x,y). 

In Section 12.11 we recall some classical definitions and propositions needed in the se- 
quel, we introduce the particle system with more details and we state the main result. 
Theorem 12.41 In Section 12.21 we derive the conditions on several examples (multitype 
contact processes, conservative dynamics and reaction-diffusion processes); we also detail 
the conditions on models with transitions of at most one particle per time. In Section 
I2.2.5l we apply the attractiveness conditions and the so-called u-criterion technique to im- 
prove an ergodicity result for a model of spread of epidemics. Others applications to the 
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construction of non-trivial invariant measures in metapopulation dynamics models will be 
presented in a subsequent paper (see f3|). In Section[3]we prove Theorem l2.4l the coupling 
is constructed explicitly through a downwards recursive formula in Section 13. 2| where a 
detailed analysis of the coupling mechanisms is presented. We have to mix births, jumps 
and deaths in a non trivial way by following a preferential direction. Section [J] is devoted 
to the proofs needed for the application to the epidemic model. Finally we propose some 
possible extensions to more general systems. 

2 Main result and applications 

2.1 Stochastic order an attractiveness 

Denote by 5 = Z'^ the set of sites and let X C N be the set of possible states on each site 
of an interacting particle system {rit)t>o on the state space Q = , with semi-group T{t) 
and infinitesimal generator C given, for a local function /, by 
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x,yeSa,l3eX k>0 

+ (<' + P'KfiS'yV) - fiv)) + (^;'/ + Pa'){fiS^'v) - /(r?))) 
where Xa^g is the indicator of configurations with values (a,/3) on (x,y), that is 



1 if r]{x) = a and r]{y) = j3 
otherwise 



and Sx^''^ : Sy, Sy^, where k > 0, are local operators performing the transformations 
whenever possible 

{r]{x) — k ii z = X and r]{x) — k & X, r]{y) + k G X 
7]{y) + k iiz = y andr]{x)-keX,7j{y) + keX (2.2) 

r]{z) otherwise; 

(S^n)(z) = I ^ 'if z = y and r]{y) + keX 

^ 1^ r]{z) otherwise; 

(S-''v){z) = I ^'^^^ ~ ^ if ^ = y and r/(y) -keX 
y \ r]{z) otherwise. 

The transition rates have the following meaning: 

- p{x, y) is a bistochastic probability distribution on 7J^; 

- ^p(a^, y) is the jump rate of k particles from x, where 'q{x) = a, to y, where 

^(y) = 

- i?|^'^p(a;, y) is the part of the birth rate of k particles in y such that ri{y) = (5 which 
depends on the value of rj in x (that is a); 

- R~^^ifp{x, y) is the part of the death rate of k particles in x such that ri{x) = a which 
depends on the value of rj in y (that is /3); 



4 



- Pg is the birth/death rate of k particles in r]{y) = f3 which depends only on the 
value of r] in y: we call it an independent birth/death rate. 

We call addition on y (subtraction from x) of k particles the birth on y (death on x) or 
jump from x to y of k particles. By convention we take births on the right subscript and 
deaths on the left one: formula ()2.ip involves births upon /3, deaths from a and a fixed 
direction, from a to f3, for jumps of particles. We define, for notational convenience 

n°'J := + P',; n;5° := i?;^/ + p-^ (2.5) 

We refer to for the classical construction in a compact state space. Since we are 
interested also in non compact cases, we assume that {r]t)t>Q is a well defined Markov 
process on a subset C VL, and for any bounded local function / on 

V?? G flo, lim ^(^)-^(^) ~ ■^(^) = Cf{r^) < oo. (2.6) 
t^o t 

We will be more precise on the induced conditions on transition rates in the examples. We 
state here only a common necessary condition on the rates. 

Hypothesis 2.1 We assume that for each {a, 13) € 

N{a, /?) := sup{n : T^^ + n°'} + n^^'^ > 0} < oo, 

E(rl/3 + n°'J + n;J'°)<oo. 

k>0 



In other words, for each a, /3 there exists a maximal number of particles involved in birth, 
death and jump rates. Notice that N(a,P) is not necessarily equal to N{f3,a), which 
involves deaths from /3, births upon a and jumps from /? to a. 

The particle system admits an invariant measure /x if /i is such that P^{r]t & A) = fJ-{A) for 
each t > 0, ^ C where is the law of the process starting from the initial distribution 
fi. An invariant measure is trivial if it is concentrated on an absorbing state when there 
exists any. The process is ergodic if there is a unique invariant measure to which the 
process converges starting from any initial distribution (see Definition 1.9]). 
Given two processes {Ct)t>o and {Ct)t>Q, a coupled process {£,t,Ct)t>o is Markovian with 
state space x Qq, and such that each marginal is a copy of the original process. 
We define a partial order on the state space: 

VC, C € f^, e < C ^ (Vx G S, e(x) < C(x)). (2.7) 

A set y C is increasing (resp. decreasing) if for all ^ G y, r/ G such that ^ < r] (resp. 
^ > 77 ) then 7/ G y. For instance if G ri, then = {r/ G ^2 : ^ < r/} is an increasing set, 
and = {77 G ^2 : ^ > r?} a decreasing one. 
A function / : ^ M is monotone if 

Ve,r?Gf^, ^<V^ f{0<fiv)- 

We denote by A4 the set of all bounded, monotone continuous functions on Cl. The partial 
order on Q induces a stochastic order on the set V of probability measures on Q endowed 
with the weak topology: 

yiy,iy'GV, iy<iy' ^{yf eM,iy{f)<u'{f)). (2.8) 

The following theorem is a key result to compare distributions of processes with different 
generators starting with different initial distributions. 



5 



Theorem 2.2 Let {Ct)t>o cmd {Ct)t>o be two processes with generators C and C and semi- 
groups T{t) and T{t) respectively. The following two statements are equivalent: 

(a) f eM and < Co implies f (t)/(^o) < r(t)/(Co) for all t > 0. 

(b) For u,u' eV,u < v' implies vT{t) < v'T{t) for all t > 0. 

The proof is a shght modification of [111 proof of Theorem II. 2. 2]. 

Definition 2.3 A process {Ct)t>o is stochastically larger than a process {^t)t>o if the equiv- 
alent conditions of Theorem \2.2\ are satisfied. In this case the process {(,t)t>o ^-5 stochasti- 
cally smaller than {Ct)t>o md the pair {£,t,Ct)t>o is stochastically ordered. 

Attractiveness is a property concerning the distribution at time t of two processes with the 
same generator which start with different initial distributions. By taking T = T, Theo- 
rem [2?2] reduces to [111 Theorem II. 2. 2] and Definition 12.31 is equivalent to the definition 
of an attractive process (see [111 Definition II.2.3]). If an attractive process in a compact 
state space starts from the larger initial configuration, it converges to an invariant measure. 

The main result gives necessary and sufficient conditions on the transition rates that 
yield stochastic order or attractiveness of the class of interacting particle systems defined 
by I^Jii. It extends jSi Theorem 2.21]. 

We denote by 5 = S{T,R,P,p) = {r^,/?, ^'J^'',p(a;, y)}{a,/3gx,fc>o,(x,y)e52}' ^he 

set of parameters that characterize the generator ()2.ip of process {f]t)t>o- We call the 
latter an S particle system and we write rjt ~ S. Given (a,/3), (7,5) € X'^ x X^, we write 
(a,/3)<(7,5) if a < 7 and /3 < 

Let S = S(T, R, P,p) and S = S{T, R, P,p) be two processes with different transition rates 
(but the same p). For all (a,/3) G X"^, let N{a,/3) be defined as N{a,f3) in Hypothesis 12.11 
using the transition rates of S. 

Theorem 2.4 Given K €z N, i := {ji}i<K, m := {'m'i}i<K, h := {/ij}j<ii:, three non- 
decreasing K-uples in N, and a, (3, 7, (5 in X such that a (3 < 6, we define 

K 

Ia:= /f(j,m)= \J{keX -.miykyd-P + ji} (2.9) 

i=l 
K 



h:= /f(j,m)= \J{keX -.j-a + miykyji} (2.10) 
1=1 
K 

Ic:= lf(h,m)= |J{fc G X : mi > /O 7-a + /ii} (2.11) 

i=l 
K 

Id:= If(h,m)= [j{k^X -.6- p + mi>k>hi} (2.12) 



i=l 



An S particle system {r]t)t>o is stochastically larger than an S particle system {(,t)t>o if 
and only if 

k£X:k>5-l3+jx kGia leX:l>ji iGlb 

E n-/ + Efl.> E <f + EA,s (2.14) 

keX:k>hi kela leX:l>^~a+hi leic 

for all choices of K < N{a, {3) V A^(7, 5), h, j, m, a < 7, /? < 5. 
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Remark 2.5 The restriction K < N{a,f3) V A^(7, (5) avoids that an infinite number of K, 
la, h, Ic, Id result in the same rate inequality. Since ^ = for each k > N{a,P), if 

K > N{a, (3) no terms are being added to the left hand side of 1^2. 13\) . and adding more 
terms on the right hand side does not give any new restrictions. A similar statement holds 
for \2.14- ), with the corresponding condition K < N^-jjd). 



Remark 2.6 IfS = S, Theorem \2.4\ states necessary and sufficient conditions for attrac- 
tiveness of S. 

We follow the approach in [8] : in order to characterize the stochastic ordering of two pro- 
cesses, first of all we find necessary conditions on the transition rates. Then we construct a 
Markovian increasing coupling, that is a coupled process (^(, Ct)t>Q which has the property 
that ^0 < Co implies 

p(«o,Co){^^ < Ct] = 1, 

for all t > 0. Here p(?"'?<') denotes the distribution of (Ct,Ci)i>o with initial state (CoiCo)- 
2.2 Applications 

We propose several applications to understand the meaning of Conditions ()2.13p - p.l4p . 
We think of (a, /3) < (7, 5) as the configuration state of two processes rjt < on a fixed 
pair of sites x and y: namely r]t{x) = a, r]t{y) = (3, ^((rr) = 7, ^t{y) = S. 

2.2.1 No multiple births, deaths or jumps 

LetlV= sup N{a,l3): 

(a,/3)eX2 

Proposition 2.7 If N = 1 then a change of at most one particle per time is allowed and 
Conditions Ii2.1cl\) and ^2.14 ) become 





+ '^a,l3 




if 13 


= 5 and j > a, 


(2.15) 








ifP 


= 6 and 7 = a. 


(2.16) 


fT-1,0 






if 7 


= a and S > P, 


(2.17) 




ff-1,0 




ifl 


= a and 6 = /3. 


(2.18) 



Proof . U P <6, then 6 - ^ + ji > 6 - ^ + ji > 1 ior all K > 0, 1 < i < K so that 1 ^ 
by definition ()2.9p . Since = 1 the left hand side of ()2.13p is null; if (3 = 6 the only case 
for which the left hand side of (I2.13P is not null is ji = 0, which gives 

fe>o keia i>o leh 

Since = 1, the value K = \ covers all possible sets la and If,, namely la = {k : mi > 
A; > 0} and Jfc = {7 — a + mi > Z > 0}. If mi > 0, we get (|2.15p . If 7 = a and mi = we 
get ()2.16p . One can prove (I2.17P in a similar way. □ 
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If 13 = 5 and 7 > a, Formula (|2.15|) expresses that the sum of the addition rates of the 
smaller process on y in state /3 must be smaller than the corresponding addition rates on 
y of the larger process on y in the same state. If /3 = 5 and 7 = a we also need that the 
birth rate of the smaller process on y is smaller than the one of the larger process, that is 
()2.16p . Conditions (|2.17p - (|2.18p have a symmetric meaning with respect to subtraction of 
particles from x. 

Remark 2.8 If S = S, when = 7 and (3 = 5 conditions \2.1b]) . 112.1 8\) are trivially 
satisfied, and we only have to check 112. 15\) when a < 7 and ^2.17^ when f3 < 6. 

Proposition 12.71 will be used in a companion paper for metapopulation models, see If 
= for all a, /3, k, the model is the reaction diffusion process studied by Chen (see 
[1]) and the attractiveness Conditions (12.151) . ()2.17p (the only ones by Remark l2.8p reduce 
to 

ri^<rl(, if7>«; r^,5>ri,^ if5>/3. 

In other words we need ^ to be non decreasing with respect to a for each fixed /3, 
and non increasing with respect to /3 for each fixed a. In [3], the author introduces 
several couplings in order to find ergodicity conditions of reaction diffusion processes. All 
these couplings are identical to the coupling H introduced in Section 13.21 (and detailed in 
Appendix 13 if = 1), on configurations where an addition or a subtraction of particles 
may break the partial order, but differ from T-L on configurations where it cannot happen. 



2.2.2 Multitype contact processes 

If F^_^ = F^_^ = 0, for all {a, /3) G X^, /fc > 0, that is when no jumps of particles are present, 
all rates are contact-type interactions. Such a process is called multitype contact process. 
Conditions (f2l3]) - (l2li]) reduce to: for all (a,/3) G X^, (7,5) G X^, (a,/3) < (7,5), hi > 0, 
Ji > 0, 

fc>(5-/3+ji l>ji k>hi l>^-a+hi 

Many different multitype contact processes have been used to study biological models. 
We propose some examples with the corresponding conditions. Since the state space 
Q = {0, 1, . . . , M}^ (where M < 00) is compact, we refer to the construction in [llj . 

Spread of tubercolosis model ([17J). Here M represents the number of individuals in 
a population at a site x G Z^. The transitions are: 

= 0/?l{O</3<Af-l}, Ra^f} = 2c?Aa]l{^=o}> 

= ]1{1</3<M}, Pix,y) = -^^xr^y}- 

where y ~ x is one of the 2d nearest neighbours of site x. 

Given two systems with parameters (A, (j), M) and (A, (j)., M), the proof of [TTl Proposition 
1] reduces to check Conditions (j2.19p : 

</'/311{0</3<A/-i} + 2dXal^fj=o} <#]1{o<5<m_i} + 2dXjl{s=o}, if (3 = 6,ji = 

]l{l<a<Af,a>hi} >l{l<T,<M,^>^-a+hi}' ff 7 > a, /li > 
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which are satisfied if A < A, < and M < M. 

In the following examples we suppose S = S, that is we consider necessary and suffi- 
cient conditions for attractiveness. 

2-type contact process ( |14j ). In this model M = 2. Since a value on a given site does 
not represent the number of particles on that site, we write the state space {A, 3,0}^"^. 
The value B represents the presence of a type-B species, C the presence of a type-C species 
and A an empty site. If ^ = 0, S = 1, C = 2 then the transitions are 

-^a,/3 = 2(iAi]l{Q,=i^^=o}> -Rq,^ = 2dA2]l{Q=2,/3=0}) 

By taking hi = 0, Condition ([2J9]l is 

(l{A:=l}2dAi]l{Q,=i_^=0} + ^{fc=2}2dA2ll{Q,=2,^=0}) 

k>S-l3 

< 2dXil{j=i^s=o} + 2(iA2]l{^=2,<5=o}; 

By taking /S = 0, 6 = 1, a = j = 2 we get 2dA2 < 0, which is not satisfied since A2 > 0. As 
already observed, see [201 Section 5.1], one can get an attractive process by changing the 
order between species: namely by taking ^ = 1,5 = and C = 2 the process is attractive. 



2.2.3 Conservative dynamics 

If nj^ = 0, for all (a,/3) £ X^,k (£ N, we get a particular case of the model introduced 
in [8], for which neither particles births nor deaths are allowed, and the particle system is 
conservative. Suppose that in this model the rate ^(y — x) has the form ijp{y — x) 
for each k, a, (3. Necessary and sufficient conditions for attractiveness are given by [U 
Theorem 2.21]: 

^ Tipp{y -x)<Y^ r^,5P(y - x) for each j > (2.20) 

k>5-f5+j l>j 

r^^^p(y -x)> ^ r^_5p(y - x) for each h > (2.21) 

k>h l>^-a+h 

while (f2J3]l - (f2H]l become 

^ Vlpp{y -x)<Y, ^\MV - ^) (2-22) 

J2 ^iMv - ^) > E ^\My - ^) (2-23) 

for all /q, If), Ic and given by Theorem 12.41 

Proposition 2.9 Conditions i fEI^) -i fEI7]) anc? Conditions are equivalent. 
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Proof . By choosing jj = j, hi = h and nii = m > N{a,/3) V N{'y,5) for each i in 
Theorem Eai Conditions (fr22]l - (l2:23]l imply l^20ij-<^21^ ■ For the opposite direction, by 
subtracting (|2.2ip with h = mi > 6 - + ji to (|2.20p with j = ji 



r'^Mv-^)^ E r>(y-x) (2.24) 

mi>k>S~P+ji ■~f-a+mi>l'>ji 

and we get ()2.22p by summing K times (12 .24^ with different values of ji, nii, 1 < i < K . 
Condition (|2.23p follows in a similar way. □ 



2.2.4 Metapopulation model with Allee effect and mass migration 

The third model investigated in [3j is a metapopulation dynamics model where migrations 
of many individuals per time are allowed to avoid the biological phenomenon of the Allee 
effect (see [1], [H]). The state space is compact and on each site x € Z"^ there is a local 
population of at most M individuals. Given M > Ma > 0, M > N > 0, cf>, cf>A, A positive 
real numbers, the transitions are 

= /?11{/3<A/-1} Pj^^ = /5('/'A%<M^} + 0]l{Af^</3}) 
_ f A a - k > M - N and P + k < M, 
"''^ \ otherwise. 

for each a,/3 E X, and p{x,y) = ^Ijy^x}- In other words each individual reproduces 
with rate 1, but dies with different rates: either 0^ if the local population size is smaller 
than Ma (Allee effect) or (j) if it is larger. When a local population has more than M — N 
individuals a migration of more than one individual per time is allowed. Such a process 
is attractive by [5i Proposition 5.1, where N and M play opposite roles], which is an 
application of Theorem 12.41 



2.2.5 Individual recovery epidemic model 

We apply Theorem [23] to get new ergodicity conditions for a model of spread of epidemics. 
The most investigated interacting particle system that models the spread of epidemics is 
the contact process, introduced by Harris [TO]. It is a spin system {'qt)t>o on {0, 1}^ ruled 
by the transitions 

r/t(x) = — )■ 1 at rate f]t{y) 

r]t{x) = 1 ^ at rate 1 
See |llj and [12] for an exhaustive analysis of this model. 

In order to understand the role of social clusters in the spread of epidemics, Schinazi [17j 
introduced a generalization of the contact process. Then, Belhadji [2j investigated some 
generalizations of this model: on each site in Z,'^ there is a cluster of M < oo individuals, 
where each individual can be healthy or infected. A cluster is infected if there is at least one 
infected individual, otherwise it is healthy. The illness moves from an infected individual 
to a healthy one with rate (p if they are in the same cluster. The infection rate between 
different clusters is different: the epidemics moves from an infected individual in a cluster 
y to an individual in a neighboring cluster x with rate A if x is healthy, and with rate if 
X is infected. 
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^v(y) = i ^"^^y^ A: = 1 and 1 < r,{y) < M - 1 (2.26) 



We focus on one of those models, the individual recovery epidemic model in a compact state 
space: each sick individual recovers after an exponential time and each cluster contains at 
most M individuals. The non-null transition rates are 

o,fc ^ f 2dXr]{x) A; = 1 and 7?(y) =0 , , 

^vix)My) \ 2dl3ri{x) = 1 and 1 < r/(y) < M - 1, ^ ^ ^ 

1 and rj{y) = 
and 1 < r]{y) 
-1 and 1 < r/(y) < M, 

p{x,y) =^%~2:} ^ach (x,y) G 5^. (2.27) 

and r^(a,) ^(j^) = for all A; € N. The rate 7 represents a positive "pure birth" of the illness: 
by setting 7 = 0, we get the epidemic model in p], where the author analyses the system 
with M < 00 and M = 00 and shows [21 Theorem 14] that different phase transitions 
occur with respect to A and (j). Moreover (O Theorem 15]), if 

AV/3<i^ (2.28) 
the disease dies out for each cluster size M. 

By using the attractiveness of the model we improve this ergodicity condition. Notice that 
a dependence on the cluster size M appears. 

Theorem 2.10 Suppose 

with (j) <1 and either i) 7 = 0, or ii) j > and /3 — A < j/{2d). 
Then the system is ergodic. 

Notice that if 7 > and /3 < A hypothesis ii) is trivially satisfied. 

If M = 1 and 7 = the process reduces to the contact process and the result is a well 
known (and already improved) ergodic result (see for instance [11^ Corollary 4.4, Chapter 
VI]); as a corollary we get the ergodicity result in the non compact case as M goes to 
infinity. 

In order to prove Theorem \2.10\ we use a technique called u- criterion. It gives sufficient 
conditions on transition rates which yield ergodicity of an attractive translation invariant 
process. It has been used by several authors (see [1], [5], [13]) for reaction-diffusion pro- 
cesses. 

First of all we observe that 
Proposition 2.11 The process is attractive for all X, (3, 7, (j), M . 

Proof . Since M = 1, then iV = 1 and Conditions (l2T3]l - (l23i]l reduce to (IZT5D . (l2T7ll . 
Namely, given two configurations € ri, ^ G with £, < r], necessary and sufficient 
conditions for attractiveness are 

if < r]{x) and Civ) = viv), 



pO,l _|_ pi ^ rO'I I pi . 

a^),ay) ^ iiy) - vix)My) ^ viyV 

P > P 

— viy)' 



If ^(y) = r^{y) > 1, 2d(3£,{x) < 2df37]{x); if C{y) = r]{y) = 0, 2dXC{x) < 2dXr]{x). In ah 
cases the condition holds since ^ < ?/. □ 
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The key point for attractiveness is that R^'^^^ ^(^y^ is increasing in r]{x). 

Given e > and {ui{e)]i(zx such that ui{e) > for ah / G X, let : X x X M+ be 
defined by 

\y-x\-l 

F,{x,y) = l{,^y} (2.30) 

3=0 

for all X, y G X. When not necessary we omit the dependence on e and we simply write 

\y-x\-l 

F{x,y) = Ijx^y} Since ^^(e) > 0, this is a metric on X and it induces in a 

j=0 

natural way a metric on Vt. Namely, for each rj and ^ in we define 

Pciri.i) := Yl nvix),ax))aix) (2.31) 

where {ct{x)}xi=x<i is a sequence such that a{x) G M, a(x) > for each x £ Z'^ and 

Y a{x) < oo. (2.32) 

Denote by 0™ := {77 G S7 : r/(x) = m for each x G Z'^} and let 7]q^ G Sl^''^ and tjq G 
(which is not an absorbing state if 7 > 0). The key idea consists in taking a "good 
sequence" {ui}i^x and in looking for conditions on the rates under which the expected 
value E(-) (with respect to a coupled measure P) of the distance between rj^^ and r/f 
converges to zero as t goes to infinity, uniformly with respect to x € S. We will use 
the generator properties and Gronwall's Lemma to prove that if there exists e > and a 
sequence {ui{e)}i^x, ui{e) > for alH G X such that the metric F satisfies 

i(£F(,?0(x),7?f (x))) < -eE(F(r/0(x),r/f (x))) (2.33) 
for each x G Z'^, then 

E(limF(r??(x),ryf(x))) =0 (2.34) 

uniformly with respect to x G 5" so that the distance /Oq,(-,-) between the larger and the 
smaller process converges to zero, and ergodicity follows. 

Condition (j2.33|) leads to 

Proposition 2.12 (u-criterion) // there exists e > and a sequence {u;(e)}jgx such 
that for any I G X 

l-i 

<piui (e) - lui^i{e) < -e n,- (e) - u{€) (A V P)2dl 
I ui{e) > 

where M(e) := maxu/(e), u_i = by convention and uq = U > 0, then the system is 
ergodic. 

Hence we are left with checking the existence of e > and positive {ui{e)}i£x which satisfy 
()2.35p . Such a choice is not unique. 
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Remark 2.13 Given U = I > 0, if ui{e) = 1 for each I € X then Condition h2. 35\) 
reduces to the existence of e > such that cpl — I < —el — (A V f3)2dl for all I € X, that is 
e<l — (/> — (AV /3)2d, thus Condition l\2.28\i . Notice that in this case 



Another possible choice is 



\y-x\-l 

F{x,y)= ^ l = \y-x 

j=0 



= U and we define (u/(e))/gx 



Definition 2.14 Civen e > and U > 0, we set uo{e) 
recursively through 



ui{e) = —(- e^Ujie) -U{Xy f3)2dl + lui^i{€)j for each I e X,l ^ 0. (2.36) 

Definition 12. 141 gives a better choice of {ui{e)}i,=Xj indeed the u-criterion is satisfied under 
the more general assumption (I2.29P . Proofs of Theorems 12.101 and 12.121 are detailed in 
Section HI 

3 Coupling construction and proof of Theorem 12.4 

In this section we prove the main result: we begin with the necessary condition, based on 
[SI Proposition 2.24]. 



3.1 Necessary condition 

Proposition 3.1 // the particle system {r]t)t>o is stochastically larger than (^t)t>o, then 
for all (q,/5),(7,5) e X^ x X^ with {a, 13) < {-f,5), for all {x,y) £ S"^, Conditions (fET^ - 
(2J4\) hold. 



Proof . Let {(,,rj) G x $7 be two configurations such that S, < rj. Let y C be an 
increasing cylinder set. If ^ € 1/ or 77 ^ 

IviO = (3.1) 

Since r/t is stochastically larger than ^4, for all t > 0, by Theorem 12.21 (or [11, Theorem 
II. 2. 2] if we are interested in attractiveness) 

{f{t)ivm < {TmvM 

since u < v' is equivalent to I'iV) < i^'{V) for all increasing sets. Combining this with 
(|XTD, 

t-'mt)ivm - MO] < t-'mt)iv){v) - mv)], 

which gives, by Assumption (|2.6|) . 

{Clvm < i^Miv)- (3.2) 
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We have, by using (|2.1|) . 

x,y£S a,l3eX k>0 

+ n°;^^(v(5S;^e) - + n-f{Ms;;'y'0 - MO)) 

x,yG5 o,/3gX fc>0 

+ n°;^^ + n;'/ • (3.3) 

We write (>Cly)(??) by using the corresponding rates of S. 

We fix y E S*, (a^, 7^, 5) G X"^ with (a^,/?) < (72,'^) for all z E S, z 7^ y, and two 
configurations (^r/) G x such that ^(z) = 02, rj{z) = 7^, for all z G S, z ^ y, 
C{y) = /3, = Thus < rj. We define the set of sites which interact with y with 
an increase of the configuration on y, 

={-^S: p{z, y) > and 5] (F^^,^ + T^^, + n°;% + H^f ,) > o}. (3.4) 

fc>0 

Denote by x = {x\, . . . ,Xd) the coordinates of each x G S. We define, for each n € N, 
Cy{n) = Cy n {z € S : Ym=i ~ yi\ — suppose Cy 7^ 0, since otherwise 

(j2.13p - (|2.14p would be trivially satisfied. Given K ^N, we fix {pl}i<K zec+ ^^^^ 
each i, z, G X and p\ < ^(z). Moreover we fix {Py}i<K such that > 5, for each i. 
For n € N, let 

= U e ^ • <(y') ^ pi Ciz) > pi, for all z G C+(n)}. 

i=l 

The union of increasing cylinder sets Iy{n) is an increasing set, to which neither ^ nor 77 
belong. We compute, using (|3.3|) . 

('C]l7^{n))(0 = E y^ E (^U£i{«.-fc>P|,/3+fc>P?;}^"-/3 + ^U£i{/3+fc>Pi,}^°'^/3) 

zeC+(n) k>0 

+ E ^^(-'y)Eiu-.{^+fc>Pu(^-''^+^S 

z&C+{n) '>0 

+ E p('^y)Y.^[jtAs+i>P'y}i^'^^^s+K's)- 

z^C'+{n) '>0 

So, by setting 

K 

J{p„a,b) := J{{pi}i<K,a,b) = : a - I > pi,b + I > p^}, 

i=l 
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and by (|3.2p . if Py := minj{p^}, we get 

zeC+(n) fceJ(pz,az,/3) fc>pi-/3 ^^C+(n) fc>pJ-/3 

< E p(^^y)i E E E p(^^y) E Kf^+^l,.) 

2GC+(n) l€Jipz,"/z,S) 1>pI-S ziC+{n) 1>pI-5 

Taking the monotone limit n — )• cxo gives 

E( E E n°'>(.,.)<E( E E 

By choosing the values = a, 7^ = 7, = p^, since YlizPi^iV) ~ 

E f + E ni^. ^ E rf,, + E (3.5) 

k£j{pc,a,l3) k>pl-l3 l€J{pc,-y,S) 1>pI-S 

A similar argument with 

C- ={zeS: p{z, y) > and E + r^,5. + + n-jf) > 0} 

subsets C~{n),n G N, {P2}2gc'-(n),j<iC ^ {Px}i<x ^ ^ withp^ < a, for each i, pi > 6z, 
for all z G Cy{n), i < K, decreasing cylinder sets 

K 

Dx{n) = IJ {C G : C(x) < pj, and C(^) < pI, for ah z G ^-(n)} 

i=l 

to the complement of which ^ and r/ belong, and the application of inequality (j3.2|) to 
^,r],Q\Dx{n) (which is an increasing set since it is the complement of an increasing one) 
leads to 

E E E E <f (3-6) 

where = maxj{|)^.}, 

J-{p^,a,h) := J-{{p\}i<K,a,b) = [j{l : a-l <pi,b + l <p\}. 

1=1 

Finally, taking Py = 6 + ji + 1, p^^ = a — nii in (13. 5p . pi. = a — hi — 1, pl^ = 6 + rrii in (13. 6p 
gives ([27[3]) - (l2l^ . □ 

3.2 Coupling construction 

The (harder) sufficient condition of Theorem l2.4l is obtained by showing (in this subsection) 
the existence of a Markovian coupling, which appears to be increasing under Conditions 
(l2l3]l - (f2Ti]) (see Subsection ES]) • Our method is inspired by P Propositions 2.25, 2.39, 
2.44], but it is much more intricate since we are dealing with jumps, births and deaths. 
Let ~ 5 and ry^ ~ 5 such that 1^ Vt- The first step consists in proving that instead 
of taking all possible sites, it is enough to consider an ordered pair of sites (x, y) and to 
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construct an increasing coupling concerning some of the rates depending on r]t{x),r]t{y) 
and £,t{x),^tiy) (remember that we choose to take births on y, deaths on x and jumps from 
X to y) and a small part of the independent rates (by this, we mean deaths from x with 
a rate depending only on rjt{x) and (,t{x) and births upon y with a rate depending only 
on r]t{y) and S,tiy))- We do not have to combine any "dependent reaction" ^'^^(^x) ^t{y) 
jumps rate , ^ ^ / ^ on w with any rate R ' , . , . or T' , ^ / s if z is different from x. 

Definition 3.2 For fixed {x, y) € S"^ , for all r] ^ Q and € N /et 



\ otherwise; ^^^^ y otherwise, 

and 



j{z,w) 



p{x, y) if z = X and w = y, 
otherwise. 



An ordered pair of sites {x, y) is an attractive pair for (5, S) if there exists an increasing 
coupling for {^,t,'nt)t>o where ~ 3(^,R,h{P),q), r]t ~ S{r,R,h{P),q). For notational 
convenience we call these new systems Sf^x,y) ^(x,y) ■ 

Notice that S(^x,y) / ^{y,x)i because we take into account births on the second site, deaths 
on the first one, and only particles' jumps from the first site to the second one. The same 
remark holds for In other words in order to see if a pair is attractive, we reduce 

ourselves to a system with only part of the rates depending on the pair, and a part of the 
independent rates depending on p{x,y) (P^^^^^p(x, y) and P~^^-^p{x , y)) . 

Proposition 3.3 The process r]t ^ S is stochastically larger than ~ 5 if all its pairs 
are attractive pairs for {S,S). 

Proof . If for each pair (x, y) we are able to construct an increasing coupling for {S(^x,y) ) ^{x,y))^ 
we define an increasing coupling for (5, S) by superposition of all these couplings for pairs. 
Indeed: it is a coupling since by Definition 13.21 the sum of all marginals gives the original 
rates; it is increasing since each coupling for a pair is increasing. □ 

From now on, we work on a fixed pair of sites (x, y) € S"^ with p{x, y) > and on two 
ordered configurations {rit,(,t) € 0,^ for a fixed t > 0, r^t ~ 5, ~ 5, < r^t and 

^t{x) = a, ^t{y) = P, mix) = 7, r]t{y) = S. (3.7) 
We denote by N := N{a, /3) V A^(7, S) and by p := p{x, y). 

Remark 3.4 With a slight abuse of notation, since rates on (a, /3) involve S and rates 
on (7, 6) involve S, we omit the superscript ~ on the lower system rates: we denote by 



The rest of this section is devoted to the construction of an increasing coupling for (x,y) 
and (a, /3) < (7, 6). 
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Definition 3.5 There is a lower attractiveness problem on (3 if there exists k such that 
j3 + k > 6 and n|^'^ + > 0; /? is fe-bad and k is a bad value (with respect to P). 
There is a higher attractiveness problem on 7 if there exists I such that j — I < a and 
n~'^° + ri^ ^ > 0; 7 is Z-bad and I is a bad value (with respect to j). Otherwise /3 is k-good 
(resp. 7 is l-good). There is an attractiveness problem on (a,/3),(7,5) if there exists at 
least one had value. 

In other words we distinguish bad situations, where an addition of particles allows lower 
states to go over upper ones (or upper ones to go under lower ones) from good ones, where 
it cannot happen. Notice that Definition 13.51 involves addition of particles upon (3 and 
subtraction of particles from 7. If we are interested in attractiveness problems coming 
from addition upon a and subtraction from 5 we refer to (/3,a), ((^,7). 

We choose to define a coupling rate that moves both processes only if we are dealing 
with an attractiveness problem, otherwise we let the two processes evolve independently 
through uncoupled rates. Conditions l\2.13\} — l \2.14\) do not involve configurations without 
an attractiveness problem, so a different construction for them does not change the result. 
Since N is finite, we can construct the coupling by a downwards recursion on the number 
of particles involved in a transition. 

Our purpose now is to describe a coupling for {S(^^^y^,S(^j.^y-^), that we denote by 'H{x.,y) 
(or simply %), which will be increasing under Conditions (I2.13|) — (I2.14|) . 

First of all we detail the construction on terms involving the larger number N of parti- 
cles and we prove that under Conditions i2.1!l^) - [2.14\ ) none of these coupling terms breaks 
the partial order: this is the claim of Proposition \3.18l 

Remark 3.6 By Hypothesis \2.1\ at least one of the terms n|^'^, n^'^, n"^'*^, n"^'*^, 
and ^ is not null. We assume all these terms (and the smaller ones) positive. 
Otherwise the construction works in a similar way with some null terms. 

Definition 3.7 Let 

N+--N-S + P; N-:=N-j + a. (3.8) 

If there is a lower attractiveness problem on f3, then (3 + N > 6 {N^ > 0) and an addition 
of N particles upon /3 breaks the partial order. Such a problem comes both from birth 
(Il^^p) and from jump (r^^p) rates. If K = 1, ji =N — 6 + (3 — 1 = — 1 and mi = N 
then Condition ()2.13p writes 

1>N+ 

Notice that if I > , then /3 + N < 6 + 1, and additions of N particles upon /3 and of / 
particles upon 6 do not break the partial order on y. The construction consists in coupling 
the terms on the left hand side (which involve particles and break the partial order) to 
the ones on the right hand side (in such a way that the final configuration preserves the 
partial order on y and on x) by following a basic coupling idea. We couple jumps on the 
lower configuration with jumps on the upper one and births on the lower configuration with 
births on the upper one. Only if this is not enough to solve the attractiveness problem, 
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we mix births with jumps. 

If there is a higher attractiveness problem on 7, then j—N < a (N^ > 0) and a subtraction 
of N particles from 7 breaks the partial order. In this case the problem comes from H^^'^p 
and ^^gP'-i we use a symmetric construction starting from Condition (I2.14p with K = 1, 
h^=N-j + a- l = N'-l and nii = N: 

k>N- 

We denote by H^'^^ ^ (resp. H'^'^'^^ ^) the coupling terms which involve jumps of k (resp. 

/) particles from x to y on the lower (resp. upper) configuration; H^'^^''^ ^ (resp. H'^'^'l^ ^) 
are the coupling terms concerning births of k (resp. /) particles on y on the lower (resp. 
upper) configuration and -^"^'^'5 (resp. H'^'^^^'^g) are the symmetric ones for death rates. 

For instance -ff^'^''^''^ combines the jump of k particles from x to y on the lower configura- 
tion and the birth of / particles on y on the upper one. 

The coupling construction takes place in three main steps. 

Step 1) Suppose both f3 and 7 are A^— bad; if this is not the case one of them is good and 
the construction works in an easier way. 

We begin with jump rates. We couple the lower configuration A^— jump rate T^^p with 
jumps on the upper configuration. We first couple it with ^p, because a < a — N < 
^ — N and the final pairs of values are {a — N^P + N) < ('y — N,S + N). Let 

ttN,N,N,N fj.N A pAf \„ f-i-i-i^ 

Then if the lower attractiveness problem is not solved, that is H^'^^^'^ = T^^p, we 
have a remainder of the lower configuration jump rate that we couple with the upper 
configuration jump rate with the largest change of particles left, — 1. We go on by 
coupling the new remainder of T^^p, if positive, with F^^ gp at step. The final pairs 
of values we reach are {a — N,f3 + N), ('j — 1,6 + I), which always preserve the partial 
order on x since a — N < 'j — I when / < N. The partial order on y is preserved only if 
/3 + N < 6 + I, that is if Z > N^. For this reason we stop the coupling between jumps at 
step A^'*': this is the meaning of formula (13. 9p . 

More precisely, when F^^ > T^g and A^ — 1 > A^"*", we get the second coupling rate 

ttN,N,N-1,N-1 ._ fj.N „ ttN,N,N,N^^ pAf-l„ _ (-nN -piV . r^-l„ 

Then either H^'^^^f^'^ ^ = (F^^ — F;^ g)p, and with the two steps I = N and Z = A^ — 1 
we have no remaining part of j^p, or H^'^^^^^'^^^ = T^^^p and we are left with a 
remainder (F^^ — ^ — T^J^ )p, to go on with I = N — 2, . . . Therefore we can define 
recursively starting from ()3.1ip 

l'>l 

where J^'^ is the remainder of the jump (hence notation J ) rate F^^p left over after the 
l^^ step of the coupling construction; (|3.12p means that we couple the remainder from the 
(/ + step, J^"^^' ^7,5^^ ^*^p- 
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We proceed this way until either we have no remainder of ^p, or we have reached 
with the remainder 

tN,N+ -pAT \ ^ ttN,N,1,1 r^N \ " -r^l ^ n io\ 

•Ja^h = - 2^ Ha,k^:5 = ^a,(SP " 2^ ^^,5P > (3.13) 

1>N+ 1>N+ 

Note that, since (3 + N < 5 + 1 for / > A^"*", none of the coupled transitions until this 
point have broken the partial ordering. Proceeding until A^^ — 1 would break the partial 
ordering. Therefore we stop at the construction in Step 1 and we will couple the 
remainder with upper birth rates at Step 3a. 

Step 1 is detailed in Tables [1] and [21 where N'^'^ corresponds to the first / (going 
downwards from N) such that the minimum in (13.120 is J^i'^^- 



Table 1: N jump rate, iV^+ > N+ {J^f^ = 0) 



I 


jN,l+l . -pi „ 


ttN,N.I,1 


N 


X 




Arrf+ + 1 


X 






X 


l'>N'i+ 


N'^+ - 1 


X 




iV+ 


X 






Table 2: N jump rate, N'^+ = iV+ - 1 {J^f^ > 0) 







ttN,N,1,1 


N 


X 










iV+ 


X 





We have to distinguish between two situations: 
• if the minimum given by H^'^^^j is always the second term, we do not reach N'^^ at 

iV+ yet. Since we have decided to stop the coupling construction at step A^"*", we need 
to couple with birth rates (Step 3a) in order to solve the attractiveness problem (Table [2] 

when J^f^ > 0) and we put N'^+ = N+ - 1. 
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• if A^'^+ > N+ then = by (I3J2D . Therefore there is no need to 

continue a couphng involving T^^p, since the attractiveness problem is solved. In this 
case H^'p^'^^ = for < I < N"^^ by definition, we do not need Step 3a and we define 

In both cases, we define 

Kf,i:s--=^ forO</<iV+. (3.15) 
since these terms could break the partial ordering. Moreover, since /3 is A^— bad, we define 

•= ° ^^^^ ^ > °- (3-^^) 

Notice that such terms are not a priori null. Indeed if 7 is Z— bad and /3 is A^— good one 
can use this coupling term to solve the attractiveness problem of 11^ ^f^p. 

If there is a higher attractiveness problem, we repeat the same construction for the 
coupling terms involving the jump rate T^^p starting from (13. IIP and we define a value A^*^" 
analogous to the previous N'^'^ . The recursive formula symmetric to (I3.12p if g > 
is 

k'>k 

where J^'^ represents the remainder of the jump rate gp left over after the k^^ step of 
the coupling construction. We need to couple the remainder of ^!^gP with the lower death 
rates in Step 3a if 

4s''' = ^> - E Hti^f = r> - E > 0. (3.18) 

k>N- k>N- 

In this case we put A^^^ = A^^ — 1. By symmetry we define H^'^'^'^ = for each 

< A: < iV", H'^^^^f^ = for each A; > if N'^' > N' and H^'^^'^f = for each A; > 
if 7 is A^— bad. 

Remark 3.8 By either A^'^+ = N or N"^^ = N, that is either the lower or the 

higher attractiveness problem given by the jump of N particles is solved by the first coupling 
term. 

If there is no lower (higher) attractiveness problem we put A^"'^ = A^ + 1 (A^"'^ = A^ + 1). 
If either /3 or 7 is A^— good, we use only one of the previous constructions. Suppose for 
instance that 7 is A^— good: then the construction involving T^^p works in the same way, 

but the symmetric one is not required and we use the coupling terms H'^'^^^g only to solve 

the lower attractiveness problems induced by either T^isP (^^ Step 1) or by nj^'^p (at Step 

36). Therefore H^^'^f' (defined at Step 3a) might be non null, but we define 

H'f^y := for each / > 0. (3.19) 
If /? is A^— good a symmetric remark holds. 



20 



Step 2) Suppose /3 is A^— bad. The birth rate H^'^p could break the partial order on /3. 
We work as in Step 1 and we begin with the coupling term 

<J°:f (3.20) 

If the attractiveness problem is not solved we couple the remainder of H^^p with the birth 

rate of the upper configuration with the largest change of particles, Il^'^~^p and going 

down we couple it with nj^'^^p at /*'*-step, until / = . We define recursively starting 
from p.20p the terms 

l'>l 

where B^'^ is the remainder of the birth (hence notation B ) rate H^'^p left over after the 
l^^ step of the coupling construction. 

We proceed as in Step 1 for transitions involving births: while the minimum in ()3.2ip is 
the second term we go on downwards in L As soon as the minimum is the first term, we 
have no remainder of n|^'^p, so the lower attractiveness problem is solved and we define 

A^'^ to be this first such Otherwise we have reached A^^ with the remainder 

1>N+ 1>N+ 

We have to distinguish between two situations: 

• if the minimum given by -f^^'^'^'^ is always the second term, we do not reach at A^'^ 

yet. Since we stop the coupling construction at step A'""'", we need to couple the remainder 
with jump rates (Step 36) in order to solve the attractiveness problem (Table H] when 

B^f^ > 0) and we put = N+ - 1; 

.ifN^> N+ then B^f = (see Table [3]) . Therefore there is no need to continue a 
coupling involving H^'^p, since the attractiveness problem is solved. In this case -ff^' ^'^5 = 
for A^+ < I < A^'^ by definition, we do not need Step 36 and we define 

^alft ■= for each / > 0, (3.23) 

In both cases, we define 

Ha!^:^i - for < / < iV+. (3.24) 
The attractiveness problem coming from either birth or death rates is solved. Indeed 
Lemma 3.9 Under Condition ^2.13\) . 



Proof . If this is not the case, by (f3T7D . Table n and definition of N'^+ 

O - E <SP = Jit^' > 4: p. (3.25) 

l'>N+ 

By (lOTT) . Table Hand definition of A^^ 

"a^p- ^^,5P = B^;p >n^;5 p- (3.26) 

l'>N+ 

The sum of (l3:25]l and ([3:26]) contradicts ([33]). □ 
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Table 3: N birth rate, > iV+ {B^f^ = 0) 



I 






_ ttO,N,0,1 


N 




X 




+ 1 




X 






X 






X 







N+ 


X 







Table 4: TV birth rate, iV-B = iV+ - 1 (S^jf ^ > 0) 



a, 13 



A ttO,( 



H' 



0,N,0,l 



iV+ 



T-rO,/ 



If /3 is TV-good, then H^^':^'^s = and we define = N'^+ = N + 1. 

Remark 3.10 // 7 is N—bad, the construction involving death rates works in a similar 
way with symmetric definitions of and of a remainder D^'^ . As in Lemma \3.9l under 
Condition \2.1^^ , 

N^\JN'^->N-. (3.27) 
7/7 is N-good, then H'^^p !^f = and we define = N'^- =N + 1. 

Step 3a) Suppose (3 is iV-bad and iV'^+ = iV+ - 1. We come back to Step 1, where even 
if the remaining part of was still positive at step N~^, we decided to stop. We refer 

to Table E 

By Lemma E31 > N+. We use the upper configuration birth rate remaining from 
Step 2 in order to solve the attractiveness problem: we couple the remainder from Step 1 
of r^^p with the remainder from Step 2 of 11^'^ p, 

jjN,Nfi,N ._ jN,N+ rjjO,iV _ jjO.JV ^0,N . 
^a,l3,^,S ■~'^a,p ^^^,5 P ^^a,l3P\- 

Then, if the minimum is the second term we proceed downwards in / with terms 

l'>l 
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where JT^ ^ is the remainder of the jump (J) rate left over after step. Notice 

that j^j*^ = J„"£^ 

We proceed with (|3.28p until the minimum is the first term, in which case there is no 
remainder of T^^p, so the attractiveness problem is solved, and we define N'^^ to be the 
first such I. In other words when / > N'^^ , the coupling term is the remainder of the upper 
configuration Z-birth rate, when I = N'^^ it is the remainder of the lower configuration 
A^-jump rate, and when / < A^'^^ the coupling terms are null. 

Remark 3.11 If I > N"^^ then the minimum of ^3. 2^) is the second term, which depends 
on Step 1/ 

-if I > N^, by Tablem AH^p = U^^'\p and by (MW, = 0; 

-^f I = by mB, alt^i = n> - B^i-^' = E.>. - <>; 

-^f I < then B^j^' A = B^j^' = and i/™ = n^'>. 

In other words even if the minimum is the second term in 'fs j s ' could be null, when 
solving the lower attractiveness problem left no remainder ofH^'^^p. It means that positive 
coupling terms begin below . 

If J^f^ = 0, that is if N'^+ >N+,we put A^^^ = + 1. We give it the same value if /3 
is A^— good: in these cases H^'^f^'^ = for each I > 0. 



Table 5: Third step, N'^+ = iV+ - 1 {J^f^ > 0) 



I 




ttN,N,0,1 


N 


X 




N'^^ + 1 


X 






X 






X 







X 





If 7 is A^-bad and A^'^+ = N+ - 1 then A^"^" = A^ by Remark EiHl and we do not need 
Step 3a for the upper jump rate. If A^'^^ = A^~ — 1 (in this case A^"'"'' = A^), the con- 
struction works in a symmetric way, by coupling the remainder of the upper configuration 

jump rate defined in a similar way with the remainder of the death rate from Step 

2. If 7 is A^— good we put A^'^^ = A^ + l: in these cases II~ ^' = for each k > 

Step 3b) Suppose /3 is A^— bad and A^'^ = A^^ — 1. We come back to Step 2, where even if 
the remainder of H^^'^p was still positive at step A^~*", we decided to stop. We refer to Table 

El By Lemma E31 A^'^+ > A^+. We cannot couple B^'^ with the upper configuration 
jump rates remaining from Step 2 with / > 7 — a, because the final states we would reach 
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are {a, /3 + N), (j — 1,6 + I), which break the partial order. Therefore we put 

:= for each l>j-a (3.29) 

Then we couple the remainder B^'^ from Step 2 of H^^p with the upper configuration 
jump rates remaining from Step 1 that do not break the partial order in a to solve the 
attractiveness problem, that is with / < 7 — a, 

^a;£i ■■= ^ KsP - JZ^' a T\^,p\ with / = 7 - a, (3.30) 

if the minimum is the second term we proceed downwards in / with with I < j — a 
ttO,n,i,i \t3N,n+ sr^ TjO,N,i',i\ . rw tN,i+i . 1 

"f-a>l'>l 

= :B^i+' A [T\^sP - Jai^' A r\^sP] (3.31) 

where B^'j^ is the remainder of the birth (B) rate n|^'^p left over after step l^^ in Step 2>b. 
Notice that = i?^;f^ 

We proceed with (I3.3ip until the minimum is the first term, in which case there is no 
remainder of Yi'^p, so the attractiveness problem is solved, and we define N^'^ to be the 
first such I. In other words when / > N^'^, the coupling term is the remainder of the upper 
configuration Z-jump rate, when / = A^^*^ it is the remainder of the lower configuration 
iV-birth rate, and when I < N^'^ the coupling terms are null. 



Remark 3.12 /// > N^'^ then the minimum in 113.31]) is the second term, which depends 
on Step 2: 

• ifl> A^^+ then J^^^ A T^^p = T^^^^p by Tableland 

Hlf,i!s = 0; (3-32) 

• ifl = N'^+ > N+ (see TableU^, then J^j-^^ A T^^^^p = J^j^^^ and 

<Ki = ^> - = E <SP - T%p; (3.33) 

• ifl< A^^+ then J^j'^^ A T^^^p = J^j'^^ = by TableUl- hence 

h1'Z'!s = A,sP ^fO<l<N''+. 

In other words if there is no remainder of T^^ ^p after Step 2, the coupling term H^'^'!y'\ 
is null even if the attractiveness problem is not solved yet. It means that positive coupling 
terms begin below N'^'^ . 

If B^f^ = 0, that is if A^^ > iV+, we put N^'^ = N + 1. We give it the same value if /3 
is A^— good: in these cases -ff^'^'^'^ = for each I > 0. We refer to Table O 

Remark 3.13 In Step 3 we do not impose any restriction as in previous steps, the pairs 
of values {a — N, (3 + N), (7, 6 + 1) or {a, 13 + A'"), {'j — 1,6 + I) breaking the partial order 
on (3 (that is such that (3 + N > 6 + 1) that we could reach a priori will be avoided by 
attractiveness (see Proposition 1 3. 1 8\) . 
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Table 6: Third step, = N+ -1 {B^f^ > 0) 







ttO,N,1,1 


7 — a 


X 






X 






X 


<f ^ > 




X 







X 





If 7 is A^— bad and N = N — 1 then the symmetric construction ahows to define the 
index A^-^'^ of the last positive coupling term and the remainder T>^'^ of the death rate 

Yi~^'^p after Step 1 and k recursions in Step 36. Remark 13.131 works in a symmetric way. 

If 7 is TV-good we put N^'^ = TV + 1 and H^^^^^'f^ = for each A; > 0. 

Remark 3.14 In the indexes we have defined, the superscript B means birth, D death, d'^ 
diffusion (jump) on the lower configuration, d~ diffusion (jump) on the upper configura- 
tion, Bd birth and diffusion, dB diffusion and birth, Dd death and diffusion, dD diffusion 
and death. 

Once we constructed the coupling rates involving F^^p, n|^'^p, ^p and Il~'^'^p we 
move to rates involving less than particles. If /3 is (A^ — l)-bad (and/or 7 is (A^ — l)-bad) 
in order to solve the lower (higher) attractiveness problems given by T^^^p and Il^^~^p 

{T^'^^p and H^^^ ^^'^p) (or their remainders from previous steps) we repeat Step 1, Step 

2 and Step 3 for such terms. We proceed this way with the remainders of T^ and H^'^p 

sP ^^'^ ^^'^fp) goiiig downwards with respect to k (/) until /? is fc-good (7 is /-good). 

We couple an upper configuration jump rate ^p with a lower configuration birth rate 
^^^sP if P ~^k > I and a < ^ — I in order to solve a lower attractiveness problem, and with 
a lower configuration death rate lif^i^p if a > j — I in order to solve a higher attractive- 
ness problem: we cannot couple the same higher jump rate both with lower birth and death 
rates. A symmetric remark holds for lower jump rates. 

Coupling rates are given below by downwards recursive relations 
Definition 3.15 Given k < N and I < N , the possible non null coupling terms are given 
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by 

•Step 1 



jk.l + l _ -pfc \ ^ rTk,k,l',V \ ^ rTk,k,-l',0 

^a,li ■— ^ a,pP Z^^a,l3,y,5 / . "-g.B.-y.S 

I'yl I'yl 
jk+l,l ._ -p/ ^ _ irk', k', 1,1 _ Sr^ TTO,k',l,l 

k'>k k'>k 

= JaT ^ tncx>^-loTP + k>5)and 

{P + k <S + l) and{a-k<j-l) 

•Step 2 (births) 

V>1 

nfc+l,i _ -pfO,i TTO,k',0.l \ ^ TTk',k',0,l 

fc'>fc fc'>fc 

= ^a;?' A ^/ 5 + ^ > ^ + 

• S'iep S (deaths) 

T^k,l+l _ -r-f-kfi \ ^ rr-kfi-l'fl \ ^ TT-kfl,l',l' 

l'>l l'>l 
T^k+1,1 ._ jT-l,0^ TT-k', 0,-1,0 

k'>k 

TT—k,0,—l,0 7-,fcJ+l . 7-,fc+li -J- ^ ; \ 7 

^a,/3:7> = ^ ^ ifa>j-l>a-k 

•step 3a) (births) 



rk,k,0,r 

l'>k-5+l3 l'>l 
Tk,k,0,l _ ^k,l+l . rofc+l,/ _ TT(}.kfl.l 



rrk,l+l _ -pfc TTk,k,l',l' _ Sr^ TTk,k,0,V 



•step 3a) (deaths) 

^k+l,l _-pl ^ TTk',k',l,l rr-k', 0,1,1 

^^,5 - ^ ■y,SP '"a,/3,7,5 '"a,/3,7,5 

A;'>i-7+Q fc'>fc 
TT-k,0,l,l rj^k,l+l rr- fc,0,-/,Oi . ^fc+l,Z -j^ ^ , 

3b) (biHhs) 

f.k,l+l _ yrO,k TTO,k,0,l' _ \^ TTO,k,l' ,1' 

I'>k-S+I3 -i-a>l'>l 

H'J:/is = A ' - ^S;^;:.] ^fP + k>5anda<^-l■ 



•Step 3b) (deaths) 

^k+l,l _ jT-l,0^ \^ TT-k', 0,-1,0 rrk', 

k'>l--r+a 5-l3>k'>k 

= [^a;?' - ^5^!^] a ^^^'^ ^/ « > 7 - ^ and /? + ^ < ^ 
• The uncoupled terms are 

J a,0,^,5 I otherwise; ' "'^'^''^ [ J^'J otherwise; ^ ' 

i) H'/.'Z = I ^ + ^ > ii) i/O'O/'' = S^'^ (3.35) 

^ a,/3,7,5 I otherwise; ' "'^'T'-^ ^'"^ ^ ^ 



V a,/3,7,5 I otherwise; ' "'^'^'^ "''^ ^ ' 
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Remark 3.16 The uncoupled terms are given by the remainders of the original rates: this 
ensures that we actually get a coupling. For instance, if 13 + k > 6, then by ^3.34^ (i), 
(MM) and [3A5\} we get 

E^TTk,k,l,l rTk,k,0,l rTk,k,—l,0\ ry-k,! \ ^ f TTk,k,l,l . rTk,k,0,l \ -pk 

l>0 1>Q 

since J!;'^ = for each I <k-5 + f3. If/3 + k<5, by (i) and 1(31^) 

E/ TTk,k,l,l . TTk,k,0,l I TTk,k,—l,0\ jk,l j_ \ ^ ( TTk,k,l,l . Tjk,k,—l,0\ -pk 
V-"a,/3,7,(5 + -"a,/3,7,(5 + -"a,/3,7,(5 > ~ '^a,l3 + Z^^-^ a,l3,-f,S + -"a,/3,7,(5 ) ~ a.^P 

since by i3. 1 5\) -ff^'^'*^''^ = for each I > 0. One gets all the marginals by summing the 
coupling terms in a similar way. 

Remark 3.17 Coupling rates involving the largest change of N particles constructed above 
can be obtained by an explicit calculation starting from Definition \3.15[ For instance 
Formula 13. 1^) corresponds to 13.15\) with k = N. Indeed if (3 is N-bad the coupling terms 

^a^ty 's '^ = by and if k = N then Y.k'>NiHa,iJ,5 + H^a^p!^l&) = 0- Formulas 

i3.r/\j , i3.21\) . i3.28\) and Ii3.31]} can be obtained in a similar way. 

The easier formulation of coupling H when = 1 is derived in Appendix |Al 

We explicitly constructed the coupling rates such that the largest change of N particles 
breaks the partial order: the following proposition proves that coupling T-L is increasing 
for such rates, under Conditions ()2.13p and ()2.14p if /3 or /and 7 are bad values. Indeed it 
states that all coupling terms that would break the order of configurations are equal to 0. 



Proposition 3.18 i) If is N-bad and I < IV~^ then under Condition h2.13\) 

ttN,N.I,1 _ ttO,N,0.1 _ ttO.N,1,1 _ ttN,N,0.1 _ p, /o q7^ 

a) If J is N-bad and k < N~ , then under Condition i2.14\ ) 

Tk,k,N,N _ rT-k,0,-N,0 _ TTk,k,-N,0 _ tt- 



TTk,k,N,N _ TT-k,0-N,0 _ TTk,k-N,0 _ TT-k,0,N,N _ /o qo\ 



Proof i) Suppose ^ < iV+. Then = H^^^"^:!;^ = by and 

• Suppose N'^+ = iV+ - 1, then N^ > N+ by Lemma EJ] and H^^f^'^^ = for each / by 
(^M)- Since H^f^:^j is null for each / < N'^^ by Tabled we prove that N'^^ > N+. 
Assume by contradiction that N"^^ < iV+. By Definition of N"^^ , ([S^H]), Tableland 
Remark 13.111 

ttN,N,0,N+ -nO,N+ TyN,N++l . /ttO-^+ n 

=^-y,S P-\N+=NB}i^a,^P- 2^ n^.S^') < ' (3-39) 

l'>NB 
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We compute J^j^^ if I > N'^^ . By (pm]) . (pmi) and Remark EH] 

rrN,l + l _yN ^ "S^ r'' "S^ rrN,N,0,l' 

^a,l3 -^a,l3P- ^ y,5P - / ,^a.B.-y.S 

l'>N+ l'>l 

l'>N+ l'>l 

- %^^>o K,r^ - «> - E • (3-40) 

Hence by (f53U]l . 

This implies that if either A^^ = iV+ or > iV+ 

+ <>> E r?,!^- E n^i^ (3-41) 

l'>N+ l'>N+ 

which contradicts (13. 91). 



• Suppose = N+ - 1, that is B^f^ > 0. Then N'^+ > N+ by Lemma [33] and 
Kf,7J = for each I by ([311]). 

Notice that Condition p.l3p with K = 1, mi = and ji = — 1 reduces to 

<>S E E <SP (3.42) 

F>Af+ 7-a>/'>Af+ 

which contradicts B^'^ > if > 7 — a. Therefore we assume < 7 — a. 

Since -f^^^'^'^ is null for each / < N^'^ by Table El we prove that N^'^ > N^. Suppose by 

contradiction that N^'^ < N+ <j-a. By Definition of A^^^*^, ^Mi), Tableland Remark 



0,N,N+,N+ _-rN+^ jN,N+ + l ^ (r^* 

N,N++1 



=r^;p-]i^^^^^,+j(r^^p- Yl KsP)<^af^'- (3-43) 

l'>N''+ 

We compute B^j'^^ if I > N^'^. By i^Mi), and Remark EH] 

V>N+ ■y-a>l'>l 
l'>N+ ■y~a>l'>l 

=<,P- E E <sP 

l'>N+ N''+A{'y-a+l)>l'>l 

-%_„>^.+>z}(ri^rp-(r^/3P- E (3.44) 
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Hence by ^M), 



a,0 



1,5 f '^{N+=N<i+ 



l'>N'i+ 



If 7 — a > N'^'^ > IV^ this is equivalent to p.4ip . hence a contradiction. If N'^^ > 7 — q 
we get 

l'>N+ ■y-a>V>N+ 

which contradicts (|3.42|) . 

Claim {ii) is proved by symmetric arguments. □ 

Remark 3.19 As a consequence, Tahles\^and\^ (and the symmetric ones) do not contain 
any coupling term breaking the partial order between configurations. 

3.3 Sufficient condition 

We complete the proof of Theorem 12.41 by 

Proposition 3.20 Under Conditions ^2.13\) - [2.14\ ), % is increasing. 

In order to prove Proposition 13.21)1 we define a new system S (in fact a new pair of systems 

S := {S,S) by Remark |3.4|) . depending on S and H, whose rates are those of S to whom 
we subtract the coupled rates of V. involving changes of N particles. 

Definition 3.21 Given S, S has the transition rates 



r 



J^j tfik = N,(3 + N>6) 

Ja,p ifk<Nor 

{k = N,^ + N <6); 



1,N 



if{l = N,j-N<a), 



^ ,5 ifl<N or 

{l = N,-f-l>a); 



O ifk<N; 



n^'i^' =^^,6 for each I; 



^a^P =Da!f}^ for e'^c/i k; 
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V. 



U-fp ifKN. 



Our plan consists in working by induction on the largest change of particles n{S) which 
causes either a lower or a higher attractiveness problem: given the particle system S, it is 
defined by 

W ■.={k : (A: > 5 - /3 or A; > 7 - a) and {T^^ + T^,, + n°'J + H"^'" > 0)} 
n(5) =n{S, a, /3, 7, 5) := | _ ^ _ (3.45) 
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Remark 3.22 If j (resp. /3) is N-good, then n{S) = sup{A; : {k > 6 — (5) and (F^ ^ + 
n°',^^ > 0)} (resp. n{S) = snp{k : {k > ^ - a) and (F^ ^ + H"^'" > 0)};. 

By Remark 13.61 n[S) = N. Let n = n[S). 

We prove that if S satisfies Conditions (ITnil - dm]) and n{S) = N, then H = HiS) 
is increasing. The induction hypothesis is: if a particle system S* satisfies Conditions 
(^U^ - U^nM and n{S*) <N-1, then n{S*) is increasing. 

We give an outhne of the proof: suppose that the induction hypothesis is satisfied. We 
defined a new system S. By Proposition 13.231 n < — 1 and by Proposition 13.241 it 
satisfies Conditions (I2.13p - ()2.14p . Therefore we can use the induction hypothesis and 
H = "^(5) is increasing. This impfies, by Proposition 13.27] that 7^(5) is increasing. 



Proposition 3.23 // either (3 or ^ (or both) are N-bad, then n < N — 1. 

Proof. If (3 is N-had, we prove that fj^^ = = 11°'^. By Definition [MB and (fCTD i) 

_ N,l _ ttN,N,0,0 

and by Definition 13.211 and (13.350 i) 

Since (3 + N > 6 + 0, by Proposition 13.181 i) . H^'p^f ~ ^a'^'^'s ~ ^ ^^^'^ done. 

— N — —N 

In a symmetric way ,5 = 11^ 5 ' = if 7 is A^-bad. 

If both (3 and 7 are iV-bad, by ^J5\i . n < N - 1. If /3 is A^-bad but 7 is A^-good, then 

n = sup{k : {k > 5 — (3) and (F^^^ + ^J^p) > 0}. Again n < A^ — 1 since F^^ and ^^a^p 
are null. The same conclusion holds if 7 is A^-bad and f3 is good. □ 

The harder part is: 

Proposition 3.24 If S satisfies Conditions i\2.13\f ~ l[2.14h then so does S. 

Proof . We prove that for all K, h, j,m, la, lb, Ic, Id in Theorem 12. 4t 

E ^^', + EK,<EKUEKs (3.46) 

fc><5-/3+ji kela l>ji leh 

E ni'° + r;,, > Y: n;f + j:r;,,. (3.47) 

k>hi kela i>7-a+/ii leic 

We prove ()3.46p . Since, by symmetry, the proof of (|3.47p is similar, we skip it. 

Remark 3.25 Let A = {a e X : a < K,ja > N^}, then for each k > 6 - /3 + 

ji > N such that ji > ja we have F^^j = by Definitions of N and \3.21\ Therefore 
^ 



Ka 



min A if A^% 
K + 1 otherwise. 



Ka-1 Ka-I 

Ia^= U {^i>k>5-P + k], I^^= \J {7-a + mi>l> ji} (3.48) 

i=l i=l 
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then condition 



k>S-P+h k€l^A l>h l^I^A 

implies \3.4-6^ and we can suppose without loss of generality 

N+ > Jk. (3.49) 

7/7 is N-bad a similar remark involving IV^ and variables N'^~ , , N^'^ and N'^^ holds 
by symmetry. 

If /3 is iV-good then Condition (I3.46P is trivially satisfied. We suppose that both [3 and 
7 are N-had. If 7 is A'^-good the proof is similar but easier, then we skip it. Hence we 
suppose iV+ > and iV" > 0. By Proposition Km TT°'^ = and by Definition KT\\ 



^a^p = for each I < N. Therefore 



j: Tt'>= j: (3.50) 

fc><5-/3+ii k>5-l3+ji 

By Proposition I3.23( p = 0. Moreover -ff^'^V^'" = for each k>6 — f3 + ji>6 — phj 
^JE)- Then by Definition [MH dXT^ and (&]) we have 

T.<^P= E (rl.p-<Xf) (3.51) 

where /„ \ {N} is the shorthand for Ial{N0a} + \ {N})l{N&ia}- 

The right hand side of (I3.46P is given by Definition 13. 2H (I3.15P and (13.150 . Notice that 
^7,(5 = by Proposition [3?23] since 7 is A^-bad; moreover -f^^'^'^'^j = if / > 7 — a by (|3.29p : 

ETT°'' r,J-\^T^ ^-Sr- f nO'' « o-O.Ar.O,; rTN,N,0,l 
^h,SP + 2^ i 1,6? - 2^ [^^j,5P - ^a,l3,-y,5 " ^0,13,1,6 

i>ji leh i>ji 

I frl ttN,N,1,1 ttO,N,1,1 \ /o r;o\ 

By Remark [3J] either A^^ > iV+ or 7^^+ > iV+. We detail the case iV"'+ > iV+, which 
contains all the technical difficulties of the proof. The other proof (for A^^^ = A^^ — 1) is 
similar. 

Since N'^+ > N+ then H'^f^^j = each / by (|3TiD . Moreover H^^fl^^s = each 
/ > 7 - a by (I3T5I1 . Since iV" > and {7 - a > / > ji} C 4 then ^ ^aA^i = 

^{NB=N+-i} E ^^°',;^'^',!5- Moreover using Tables [3] and H 

'y—a>l>j\ 
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The first term on the right hand side (|3.52|) that we consider is 

i^^,5P-H^/,^:s)-\NB=N+-i} 2^ ^a,p:i> (3-54) 

• li > , since N+ > Jk > ji by ([Ti^ . we have > ji: by using Tableland 
(lOTD . (1331) is equal to 

E n> - B-f-' = - n°'>. (3.55) 

NB>l>j^ l>jl 

since in proof of Proposition 13.181 we obtained that A^-^*^ > A^+, again 
by (1331, we have N^"^ > jy. by Table El (CTD and (IT^ . (ICTD is equal to 

N+-l>l>jl ■y-a>l>NBd ;>7V+-1 

= ^n;'>-n°'>. (3.56) 

that is ()3.54p has the same value in both cases. The second sum on the right hand side of 
(|332|) is 

E (r>-<£;^)= E ^\,sP- E r>]i|.>^..} 

l&h\{N} i&h\{N} ieit\{N} 

- E i^ZfiP- E r^J,5p)]l{i=Ard+| 

= E ^\,sP- E r7,5^'^{;>^''+} 

l€h\{N} leIt\{N} 

E ^7,5^')^{^'*+eM{Af}}- (3-57) 

Therefore using (I333I1 . (I335I1 . (I336D and (f3371) 

E^7i^'+E^7,5^'=En°;!5^'-n°',?p+ E E r!,,5P%>^.+} 

/>ji i&h i>h i&h\{N} lahMN} 

-{^%P- E ^l,5P)\N'^+&h\{N}}- (3-^^) 

We use (133(1 . (|33B . (f338]l and Condition (|2l3]l to check that Condition (fOH]) is satis- 
fied. 

Case A: Suppose N'^~ = N. In this case T^^ > T^^^ (by (IXTTl) ) and H^^^^^f = for 
each k <N. Therefore by ([Xn]) 

Erl/3= E (3-59) 

• Suppose N'^+ ^ h. Since iV"'+ > iV+ > j/< by (fS^lQl) . we must have A^'^+ > 7 - a + m^, 
hence 4 \ {{iV} U {iV > / > N'^+}} = h- It implies that Eieh\{N} ^\,5hi>N^+} = 

by (l33nD . (I33TD and (f338]) Condition (jHIIHIl becomes 

E E rl,<En% + Er;,. 

A;><5-^+ii kGla\{N} l>ji leh 
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which holds by Condition (|2.13|) . 

• If N'^+ elb\ {N}, then 

-f - a + rriK > N'^^ > N+ > Jk, (3.60) 
hence (see Tabled]) by (l330|) . (I33T|) and (1338)1 Condition I^SMh is 

fc>5-/3+ii fcG/a\{A^} />ii ieh\{N} ieib\{N} 

= En%+ E ^\,s-rZ,- (3.61) 

Notice that 

{h \ {N}} U {/ > iV^+} = If, U {/ > iV^+}. 

• If iV e /a, (pTHTT) becomes 

E <^, + Eri,<En;i+ E 

k>S-/3+ji kaia i>ji iehu{i>N''+} 

which holds by (ITTH]) . 

• If iV ^ /a, (13311) writes 

E E rl,<En;4+ E ^'.,s (3.62) 

fc>5-/3+ii fcG/aU{Af} l>ji l£lbU{l>N'i+} 

Denote by 

Ta =Ia^ {N > I > 5 - P + {N - 5 + P)} = laU {N > I > N} 

h =IbU {N + -f - a > I > N - 6 + 13} 
then by Condition (j2.13p applied to la, h 

E nS + Erl.<En% + ErU (3-63) 

k>S-l3+ji Z>ii i^j^ 

By (13301) . lb = Ife U {iV + 7 - Q > / > iV^+j and ([333]) implies (1332]) . 



Remark 3.26 This is the key passage where a Definition of la, . . . Id as a single set instead 
of a union of several sets does not work. 

Case B: Suppose = N'^^. In this case H^'^^^'l = for each I < N, the left hand side 

of n3.57p is equal to gp. Formula (I3.5ip still holds but terms Ha^'^'s^ 

null. One works as in Case A, when A^'^^ = N, by using the symmetric construction of 

k k N N 

5 ^'^d by checking the condition in different cases. □ 
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Proposition 3.27 If H := 'H{S) is increasing, then % = T~i{S) is increasing. 

Proof . We show that all coupling rates of % breaking the partial order between configu- 
rations are null. By Proposition 13.181 the ones involving particles are null and if both 
(3 and 7 are {N — l)-good then T-L is increasing. Therefore we suppose that either /3 or 7 
is [N - l)-bad. 

We prove that 1-i and T-i differ only on rates where a change of N particles causes an 
attractiveness problem. By Proposition 13.231 such terms are null for %. 

Step 1: by Definition [3l5] ([315]) on S and Definition [MH if iV > A; > A; - (5 + /3, we get 

■Yjk,k,l,l rj^k \ ^ r-fjk,k,l' ,1' ^fc.fc,— /',0\ ■, rj^l \ ^ /^k' ,k' ,1,1 ■YjO,k' ,1,1 

^a,l3,y,5 =1^ a,l3P ~ a,/3',7,(5 + J^a,f5,j,5 )\ '^[^ j,5P ~ a,/3,7,(5 + -"a,^,7,5jj 

l'>l k'>k 

\jk,N \r^,-rjk,k,l',l' -Yjk,k,-l',0.. r N,l \ ^ ,^k',k',l,l -rjO, k' ,1,1 (:a\ 

V>1 k'>k 

Note that since /'^'^ - r^- R^'^^^^^ ^ R^^^'-^fi ^r,A jN,l _ yrO,l rrN,N,l,l _ rrO,N,l,l 

iNOte tnat SmceJ^^g - ^ a,l3P ^a,l3,f,5 +^0,13,^,6 '^"Cl -'7,5 - -^S.^^ ci,p,'y,S -"a,/3,7,<5' 

by definition (fXT5|) when Q <k,l <N -I, 

TTk,k,l,l _r jk,N Sr',TTk,k,l',l' TTk,k,-l',Ox-, . ^ jN,l ST- f TTk',k',l,l , TTO,k',l,U-, {o p,^\ 
i'>Z k'>k 

By ([OTD and ([MH]), i^JJ/^'^ and ^H^^;^';;'^ have the same recursive definition when < 
/t < TV - 1, < / < iV - 1. ' ' 

Step 2: If /3 is /c-bad (otherwise H^J^if^^^^ = -f^^'^''^''^ = for each / > 0), by Definition 13.151 
(f3T5]) and Definition [321] 

J^a,fS,'y,5 —[^'-a,l3P ~ a,/3,7,<5j [-'-^..SP ~ l-"a,/3,7,5 + I^a,l3,-f,5 )\- 

l'>l k'>k 
rrrO.fc TrO.fe.O,/' -, AT,; ,-^.0,fc',0,« ^k',k',0,L-. aa\ 

= \^a,pP - ^a,/3,7,5] ^ [^7,5 " (^«,/3,7,5 + ^«,/3,7,5 )]• f-^'^^) 
l'>l k'>k 

As in previous case, since B^j = n°''^p - - H^'^J^^J, by definition p.l5p, 

Q ^ ^ 

and H^^J^^g have the same recursive definition when 0<A;<A^ — 1, 0</<A^. Notice 

— 0,N 

that in this case they coincide even if / = A^, since we may have H^'^ p > 0. By similar 
arguments and (j3.15p we prove the symmetric result with respect to death rates. 

Step 3a: Suppose that (3 is /c-bad and J^j^ > 0, that is Step 1 was not enough to solve 
the attractiveness problem induced by for S. By Definitions 13.151 ()3.15p - (l3.15p and 

■Yjk,k,0,l _rj^k \ - -Yjk,k,V,l' •sr^^k,k,Q,l' . sr^ -fjO,k' ,0,1 sr^ -fjk' ,k' ,0,L. 

I^a,lS,y,5 —[^ a,pP ~ 2^ I^a,l3,j,5 ~ ''^ a,l3,-y,s\ ^ [^^j,5P ~ -"a,/3,7,(5 ~ -"a,/3,7,<5 )\ 

l'>k-5+l3 l'>l k'>k k'>k 

_r Tk,N \ ^ -Yrk,k,V,V V^Tr^.^.O.''! . rnA^,' ST^ -nQ,k' ,0,1 ^ -rjk' ,k' 

~\-'^OL,P ~ ^a,l3,^,S - Z^^a,l3,^,5\ /\ [^7,5 " 2^ ^a,l3,'r,S " 2^ ^ a,l3,^,5 )\- 

I'>k-S+I3 l'>l k'>k k'>k 

(3.67) 

^\r,n^ t'''^ - rk Lrk,k,N,N rTk,k,0,N . r^N,l _ yrO,l ^ rrO,N,0,l rrN,N,0,l -^^ 

-^ais'-y's ^a,i3,y,s ^ave the same recursive definition when 0</i;<A^ — 1,0<Z<A^. 



34 



The equality of recursive formulas corresponding to Step 3b is proved in a similar way. 



Since all terms are defined by the same downwards induction formula, we just need 
3 check that they coincide for the initial coupling rates. The one involving T^^^p and 
^J^p is given by (fS^M]) with I = N - 1, that is 



-YyN-1,N~1,N~1,N-1 _ jN-l,N . jN,N-l _ rpAT-l ttN-1,N-1,N,N TTN-l,N-l.-N,0-i 

. mAf-l _ ttN,N,N-1,N-1 _ ttO,N,N-1,N~1^ _ ttN-1,N-1,N-1,N-1 

and we are done. The one involving n|^'^~^p is given by (I3.66P with / = A^, that is 

^0,Af-l,0,Af T-rO.Af-l A T3^,N rrO.Af-l a FrrO.Af ttO,N,0,N ttN,N,0,N-i 
_ttO,N-1,0,N 

We prove that H^l^'^'s ^^'"^ ~ ^^a^'-y's ^^'^ death rates by symmetric arguments. 
Therefore Ti and Ti are identical for all coupling rates involving a change of less than N 
particles which cause an attractiveness problem: the conclusion includes uncoupled rates. 
The claim follows since Ti is increasing and by Proposition 13. 181 □ 



3.3.1 Proof of Proposition [37201 

We do an induction in two steps. First of all we suppose that /3 is A^-bad and 7 is A^-good, 
that is (3 + N > 5 and 'j — N > a. We skip the similar symmetric case. 
Let <S^ := 5 and suppose that Proposition 13.201 holds for each system S* such that 
n{S*) < N — 1. Notice that since 7 is A^-good, then 7 is n(5*)-good for each S* and, 
by Remark I3.22| definition of n{S*) involves for each S* only rates that cause a lower 
attractiveness problem. 

Remember that Definition 13.211 of S depends on the original system, that is 5 = S{S). 
We define S* = S{S*^^) and nj = n{S*) for j G N,i < A^ - 1. 

If f3 is (A^ — l)-bad, by Proposition I3.23"l n^v-i < A^ — 1, and 5jv-i satisfies Conditions 
()2.13p ~ ()2.14p by Proposition 13.181 We define a coupling Hn-i for S^j^-^ as in Definition 
I3.15[ and by induction hypothesis it is increasing. By Proposition 13.271 if T-Ln-i is increas- 
ing then so is Hn- 

We have to check the induction basis: we proceed downwards with definitions of the new 
systems S* until /3 + j < 5. If /? < (5, then S* with j = 5 — /? is attractive, since /? + j < 5 
and there are no attractiveness problems. If /3 = 5, then the attractiveness of a system 
with ni = 1 under Condition (|2.13p is proved in Appendix [XI 

We proved that all systems such that there is only a lower (or only a higher) attrac- 
tiveness problem are attractive under Conditions ()2.13p - (l2.14p . Now we prove that this 
holds also for a system where both (5 and 7 are A'^-bad. 

li j3 + N > 5 and j — N < a, the definition of n{S) is given by (|3.45p but everything 
works in a similar way: we define = S, S* = S{Sj^i) and rij = n'{S*) for j < A^ — 1. 
If both /3 and 7 are (A'' — l)-bad, by Proposition 13.231 n^v-i < A^ — 1 and satisfies 
Conditions ()2.13p ~ (l2.14|) by Proposition 13.181 We define a coupling Hn-i for as in 

Definition 13.151 by induction hypothesis it is increasing. By Proposition 13.271 if Hn-i is 
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increasing then so is "Hat. 

We have to check the induction basis: we proceed downwards with the definition of a new 
system S* until either /? or 7 are j-good, that is until j = ((^ — /?) V (7 — a). 
If 5 — /3 ^ J — a, then only one attractiveness problem is present and the claim follows 
from the first part of the proof. 

Ifj = (5 — /3 = 7 — a>0, then Sj is attractive since there are no attractiveness problems; 
if J =6 — 13 = 7 — Q! = 0, then the attractiveness of a system with ni = 1, /3 = 5 and 7 = 
under Conditions (I2.13p - (l2.14p is proved in Appendix O □ 

3.4 Sufficient conditions on more general systems 

In order to show the sufficient conditions of Theorem 12.41 we restricted ourselves to tran- 
sition rates on a given pair of sites {x,y): if Conditions ()2.13p ~ (l2.14p are satisfied we can 
construct an increasing coupling for the system ^.nd the final increasing coupling is 

given by superposition of couplings for all pairs of sites (see Section 13.21 and Proposition 
13. 3p . We use neither the translation invariance of P^^^y nor the fact that the smaller 
and the larger systems share the same p{x,y). Therefore we can state the result for a 
more general pair of systems G = {R'^^'^^^^f^y^{x,y),T^^^^^^(^y^{x,y), P^^^^ and 

Corollary 3.28 A particle system ijt ^ G is stochastically larger than ^ G if for each 
(x, y) € S"^, (a, /3) < (7, 5), (a, /3) G X^, (7, 6) G X"^ , with p = p{x, y) and p = p{x, y) 

k>&-p+ji keia i>ji leh 

k>hx keld l>-y-a+hi iGlc 

for all choices of K < N{a,j3) V N{'y,5), k, j, m in Theorem \2.4\ 

Such conditions are not necessary. We use Corollarv 13.281 with the comparison technique 
with oriented percolation in [3J to prove survival of species in metapopulation models. 

4 Proof of Theorem [230] 

Since F{x,y) = F{0,y — x) = F{y — x,0), with a slight abuse of notation we write F{y — x) 
instead of F{x, y). 

Proof of Proposition \2.12[ 

We treat separately the cases 7 = and 7 > 0. 

i) If 7 = then the Dirac measure 60 is invariant. We denote by P(-) the independent 
coupling measure and by E(-) its expected value. We fix x G S and we compute the 
generator on F(r/*^(x)) in Definition (j2.30p . where {ui{e)}i^x satisfies Hypothesis (j2.35p . 
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Let t > 0, we denote by rj^^ (x) = I. By (IT^ 
/:F(ryf (x)) =CF{1) = ll{,=o} [a 7?f (y)(F(l) - F(0)) 

+ ]i{i<KM-i} [w E ^t'^(y) + ^'^)(^(^ + 1) - ^(^)) + ^(^(^ - 1) - ^(^)) 



+ 1 



{«=M} 



M{F{M - 1) - F(M)) 



i{/=0} 



+ i{i<«<M-i} ^ vi' iv) + i^hi - ly'l-i] + 1 



{l=M} 



MUM- 



Ul - lui_i 



y^x 



<1<M} 



Ul - lui^i 



yr^x 



y~x 



<(A V/3)n^r/f (y) -e^Uj - {\y P)2dul. 

j=0 



(4.1) 



yr^x 



Let 7]^^ G By translation invariance and Definition (|2.30p 



|E(F(,?f(:r)))=E(/:F(7?f(x)))<E((AV/3)u5]r?f(y)-6 ^ - (A V /3)2d«77f 

= -6E(F(r,f(x))). 
By Gronwall's Lemma we get the result. 

a) If 7 > then 5o is not any more an invariant measure. Let {(,t,'nt)t>o be a coupled 
process through the basic coupling probability measure P such that ^ ^-nd tjq G . 
We fix X S Z'^. We denote by A; = ^t(x), / = Ct{x) ■= rit{x) — S,t{x) and k + 1 = rjt{x). First 
of all we prove 



(x-)-l 



Lemma 4.1 



£F(/c, /c + <A Y^^miv) - it[y))ui +l{(t)ui - ui^i) . 



(4.2) 



yr^x 
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Proof . By Definition 

^F{Ct{x),r]t{x)) =l{k=i=o} [a ^ Ct{y)ui 



+ 1 



{fe=0,0<«<M-l} 



yr^x 



+ \k=o,i=M} [{\Y,it{y)+i){-ui.i) + M{-UM.i) 

y^x 



+ 1 



y^x 



(4.3) 



We prove that 

CF{k, k + l) <%=o,o<i<A/-i} - (A X] + 7)^ii-i - ^ui-i + (/^ X] + 



yr^x 



+ 1 



{{fc=0,0<i<M-l}':} 



-1 



(44) 



y~x 



Notice that if /c = then r]t{x) = k + I = I. 

If > and k + l = M, then the last term in the right hand side of (|4.3p is smaher or 
equal to 

yr^x 

The same inequality holds when / = and k = M, that is ^t(x) = rit{x) = M, since the 
last term in the right hand side of ()4.3p is null and l(f)ui — lui^i = 0. If A; = and I = 0, 



y~x 



Therefore (g^D holds. 

Since ui is non increasing in /, 



y^x 



yr^x 



If /? < A, 



and by (j4.4p the claim follows; if /3 > A 



y^x 



y~x 



l(t)ui - lui^i + XYCt{y)ui + U/3 - A) ^?7t(y) - jjui. 



(4.5) 



yr^x 



y^x 



Since /? - A < 7/(2dM), then (/3 - A) r?t(y) - 7 < and the claim follows from ([131) 



and (1431). 



□ 
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By Lemma |4.H (|2.35|) and translation invariance 

E(£F(fc, k + l))< - lui^i + X^Ctiy)ui 



<e(^ - e^Uj - XU2dl + X'^Ct{y)U^ = -eE(^F(A:, A: + /)) . 
j=0 y^x 

and the claim follows by Gronwall's Lemma. □ 

Now we prove that Definition 12.141 satisfies the hypothesis of the u-criterion under 
()2.29p . We begin with a technical proposition. 

Proposition 4.2 Let {ui{e))i^x be given by Definition \2.14\ If (p < 1 and 

< A V /3 < 



2d 2d{l - 

then there exists e > such that ui{e) is positive, decreasing in I for each I G X and in e 
for each < e < e. 

We prove Proposition 14.21 by induction on /. We need the following lemma to use the 
induction hypothesis: 

Lemma 4.3 If for I G X, 

AV/3< 



2d{l - (pi) 

and there exists e such that Ufc(e) > for each k<l — 1, 0<e<e; then there exists 



< e* < e such that 



U< ^'-^^^ V (4.6) 
2d{X y 13) ^ ' 



Proof . We prove by a downwards induction on < /c < / — 1 that: 
if there exists < e* < e such that 

J{k, e*) := -(A V P)2dU{l + + ... + </)'=) + ui_k~i{(*) > (4.7) 

then J{k - l,e*) > 0. 

Indeed, if A: = / - 1 then is the assumption (A V /3)2(i(l - </>') < 1 - 0. 

Suppose there exists < e* < e such that J{k,e*) > 0. By Definition 12.141 (jiTTl) is 

equivalent to 



-e 



EfJ, + (i - fc - 1)( - f/(A V fl)2d + m-t-2(e')) 



that is 

l-k-2 

-(A V ^)2dU{l - A: - 1) [(1 + ... + (j)'')(l) + l]-e* ^ Uj{e*) + {I - k - l)n,_fe_2(e*) > 0. 

j=0 

Therefore 

J{k - 1, 6*) > -(A V p)2dU{l + <P+... + (/>'+^) + ui.k-2{e*) > ^ _ ^ 
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Since by hypothesis Uj{e) > for each < e < e, then J{k — l,e*) > and by induction 
holds for each < k < I — 1. By taking k = we get 



-(A V/3)2dC/ + ^//_i(e*) > 
which is the claim. □ 



Proof of Proposition \4-^ 

We prove by induction on / G X that there exists e/ such that for each < e < e/ and for 



each < 7 < /, Ujie) is positive, decreasing in j, U < — ^ ^ ^, and 



> ±u,{e) > -CuU) 



with Cu{j) the solution of 



c„a) = ^Mi^. ,4..) 



This gives 



CuU) = ^(1 + ^ + + . . . + ^-^) = U^^^ (4.10) 
hence Uj(e) > is also decreasing in e. 

We prove the induction basis when 1 = 1. By Definition 12.141 

-euo-U(Xy f3)2d + uo -e - (A V /3)2d + 1 
^^i(e = = UQ . 

Since (A V /3)2d > 1 — we can take e < ei small enough to have 1 — (A V /3)2d — e < cj), 
that is ui{e) < uq = U; since (A V /3)2d < 1 then U < U/{2d{X V /3)) and, by taking 
e < small enough, ui{e) is positive; moreover notice that ui{e) is decreasing in e, and 
^ui(e) = — ^ = —Cu{l) for each e. Hence the induction basis as well as the hypothesis 
of Lemma 14.31 are satisfied for / = 1: there exists ei = ei A such that if e < ei then 
U > ui{e) > 0, ui(e) is decreasing in e, and (14. 6p holds. 

Suppose there exists e;_i > such that for each < e < e/_i, then > ■^Uj{e) > —Cu{j), 

Uj(e) is decreasing in j for each j < I — I, U < — ^ ^, and Uj(e) > for each j < I — I. 
3\ J t> J J - 2c?(A v/3) ' 

First of all we prove that there exists e/ > such that if < e < e; then ui[e) < u;_i(e). 
By Definition [TTil 

By induction hypothesis, if e < e;_i then U > Uj{e) > for each j < I — 1 and we get 

-eEj=o-"j(g) - ^(A V /3)2(i/ + /^f-i(e) -[/(A V /3)2d/ + ?n;_i(e) 

(t)lui-i{e) (t)lui_i{e) 

1 (AV/3)2d 
<- - -^^ — — < 1 
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that is ui{e) < tf;_i(e) by (|4.1ip and we set ei = ei_i. 

By ()4.10p . Cu{j) is always positive and increasing in j. We prove that ()4.8p holds for 
j = 1. By Definition [m 

We begin with the right inequality involving the derivative of ui{e) in (14. 8p . By induction 
hypothesis, if < e < q_i by ([O]) . (|4.10p and (|4.12p on the one hand 

4)l^ui{e) >-lU- lCu{l - 1) = -Cu{l)l(t). 



and on the other hand 

j=0 i=0 i=0 ^' j=0 

Since < e < q_i and Uj{e) is positive and decreasing in e, as e approaches the first sum 
on the right hand side goes to 0, while the second sum is positive and increasing: therefore 
for each j < / — 1, we can take ei small enough so that ■^ui{e) < for each < e < e;. 
Now we prove that there exists > such that ui{e) > for each < e < e^. By 
Definition Em ui{e) > if 

i-i 

-U{\y P)2dl + lui_i{e) > e^Uj{e). (4.13) 

j=0 

By the induction hypothesis, assumptions of Lemma 14.31 are satisfied, hence there exists 
< e* < Q_i such that is satisfied. Thus -C/(A V /3)2dl + lui-i{e*) > and we can 
choose < e* such that 

-U{\ V f3)2dl + lui_i{e*) > (*ilU (4.14) 

Ife<er(<ei-i) 

-U{\ V p)2dl + lui_i{e) > - [/(A V fi)2dl + lui^i{e*) 

i-i i-i 

j=0 j=0 

which is ()4.13p . By taking < e < eg A . . . e;_i A q A , ui{e) > and is decreasing in / 
and e, and the claim follows. □ 

Proof of Theorem \2.1(A 

Let {ui{e)}i(zx be given by Definition (|2.36p . By Proposition 14.21 we can choose < e < e 
such that ui{e) > for each / G X and ergodicity follows from Proposition 12 . 1 2l □ 
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A Appendix: Explicit coupling construction \i N = 1 

We detail the coupling when = 1 to understand the simplest construction. Note that 
we have to mix birth with jump rates also in this case. 

There is a lower attractiveness problem (/3 + 1 > 5) only if /3 = (5, and a higher one 
(7 — 1 < a) only if 7 = a. Definition 13.151 becomes 

Definition A.l The non zero coupling rates ofH are given by: 
= A T\,)p ^fP = 5or^ = a■ 



Hl'XA = [n°'> - (nS A u'^j)p] A [ri,,p - (r^,, a ri,,)p] ^/ 7 > « J 



0-1,0,-1,0 _ (fr-ifi . n-i'Ou ^ 

H-^'^s' = [nS°p -iKf A n~i/)p]+ A [pi^^p - {Tl^ A ri,,)p] z/ « = 7; 

Hli-y = K,^P - K,, A Pi,,)p] A [U-}'% - (U-y A n-;/)p] ^fP<6 j 



^0,1,0,0 _ ^=^0,1 _ rrO,l,0,l _ rj-0,1,1,1 0-1,1,0,0 _ pi _ 0-1,1,1,1 _ 0-1,1,-1,0 _ 0-1,1,0,1 

a,/3,7,(5 ~ '■^a,l3P ^a,l3,'r,5 ^a,l3,-f,5^ ^a,l3,'y,S ~ ^ a,fiP ^a,l3,-y,S ^a,l3,"/,S ^a,l3,-f,5^ 

rO,0,0,l _ ■p-0,1 _ 0-0,1,0,1 _ 0-1,1,0,1 oO,0,l,l _ -pi _ 0-1,1,1,1 _ 0-0,1,1,1 _ ij-1,0,1,1 

a,/3,7,6 ~ ^'■'r,sP ^a, 13,^,5 ^a, 13,^,5^ -"a,/3,7,(5 ~ J,5P ^a,p,^,5 ^a, 13,^,5 ^a, 13,^,5 ' 

r-1,0,0,0 _ -pj— 1,0 _ 0—1,0,-1,0 _ 0—1,0,1,1 0-0,0,-1,0 _ -TT- 1,0 _ o— 1,0,-1,0 _ 0-1,1,-1,0 

a,/3,7,(5 ~ ^^a,l3 P ^a,p,^,5 ^a,l3,-/,5 ' ^a,l3,-/,5 ~ ^^-f,5 P -"a,/3,7,<5 -"«,/3,7,<5 • 
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Proposition A. 2 Definition \A.l\ gives an increasing coupling. 



Proof . The definition of uncoupled terms ensures this is a coupling. We have to prove 
attractiveness. 

Suppose /? = (5 and 7 > a. Therefore H'^/^^^^ = = ^a,^',^,^/ = 0- 

• Suppose f > r^^^. By Condition (|2.15p we must have n°'^^^p < n°'_^p, then H'^^^^^^s = 
0. Moreover 

Tj-l, 1,1,1 _ pi Tj-0, 1,0,1 _ fj-0,1 tjO, 1,1.1 _ 

by Condition ^M- Therefore = T^^^p - - fl f,p + = 0. 

• Suppose fl^p < r^ ,5. Then = f^^^p and we get = = 0. We 
have to prove that -f^^'Jj''^*^^ = 0. 

If 7 = a, by Condition (|2.16|) then -ff^'^'^j^ 5 = n[^'^p — n|^'^p A ^^^'\p = and we are done. 
Suppose 7 > a: 

[7-0,1,0,0 -pfO,! f^O,l . 1^0,1 rfi-0,1 -p;0,l . rrO,! i a f-pl pi 1 

By using the relation 

a A (c - c A 6) = (a + 6) A c - 6 A c for a, 6, c > 
with a = rl^^^p - r^^^p, 6 = n^'^p and c = n^'^^jP, we get by Condition (|2.15p . 

=n°'>-n°'> = o. 

If a = 7 and 6 > (3 we check that the coupling is increasing in the same way. □ 
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